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Abstract 

The purpose of this note is to show that the Jacobi problem of geodesies 
on ellipsoid [1], [2], [3] may be reduced to a more simple case, namely, 
to the Clebsch problem [4] . The last one is the problem with quadratic 
nonlinearity and may be solved directly by using of Weber's approach 
[5] in terms of multi-dimensional theta functions. 

1. The problem of geodesies on ellipsoid is a classical one. For two-dimensional 
ellipsoid, its solution was announced by Jacobi on 28 December, 1838 [1] (see 
Appendix). Using the remarkable substitution, he reduced this problem to 
quadratures. These results were published in the paper [2] and then consid- 
ered in details with the use of elliptic coordinates in his lectures at Konigsberg 
University in 1842/43. They can be found in the book [3]. 

Later Weierstrass described another approach and succeeded in integrat- 
ing of the equations for geodesies explicitely in terms of two-dimensional 
theta functions [6] . 

Another solution of this problem was given by Moser ([7], [8]) who in- 
cluded it to the modern scheme of isospectral deformations and described 
the relation of this problem with the geometry of quadrics and algebraic ge- 
ometry of spectral curve. For n-dimensional case, the algebro-geometrical 
approach was given by Knorer [9]. 
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Here we show that this problem may be considered without using elliptic 
coordinates and algebraic geometry as a projection of a simpler system (for 
a general description of this method, see [10], [11]). Namely, it may be 
considered as a projection of the system with a larger number of degrees of 
freedom and with the simplest (quadratic) nonlinearity, the so-called Clebsch 
system [4]. The equations of motion for this system were integrated in terms 
of multi-dimensional theta functions by Weber [5] (for modern exposition, 
see [12]). 

2. The (n — 1) -dimensional ellipsoid in n-dimensional Euclidean space R n is 
defined by the equation 

n 2 

= E v = 1 (!) 

The standard metric in Euclidean space induces the metric on ellipsoid. The 
free motion of a point particle of unit mass on such an ellipsoid is described 
by the Hamiltonian 

n = \£v!. (2) 

3=1 

The equations of motion have the form 

i 3 = yj> i/j = = -v(x,y) — , j = l,...,n. (3) 

To obtain v(x, y), we shall differentiate the equation (1) twice with respect 
to time. As a result, we get 

Y^± = (4) 

J= i j 



and 



where 



n 2 n 2 

J=l J J=l J 
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Theorem. Equations (3), (5), (6) for geodesies on ellipsoid are equivalent 
to the equations of motion for the Clebsch system, 

k 

(7) 

hk = y^(^jm ^mk — Uj rn l mk ) — (cij 1 — a fe X ) yj y k . 

m 

Here we introduced additional variables related to angular momentum tensor 

Ijk = XjVk - XkVj and u jk = — l jk (8) 

and r is a /oca/ time defined by the formula 

dr = A- 1 dt, ^ = A^. (9) 
dr dt w 

Proof. After the change of a standard time t by a local time r, the right-hand 
side of equations (3) takes a polynomial form 

-f •<•, .U-r)//,, -f//, Bfc,)^. (10) 
ar ar aj 

One can see that there are following identities: 

/c j<k j k 

(ii) 

(aj 1 - a k x ) XjX fe = ^(/ jm a; mfc - w im / mfc ) - (aj 1 - aj 1 ) yjy k . 

n 

From here, it follows equations (7). They are the equations with the simplest, 
namely, quadratic nonlinearity and they are exactly equations of motion for 
the Clebsch system [4]. 

3. It is easy to check that equations (7) have n independent integrals of 
motion, 



1 

aj — a k 



Fj = Vj + - j,k = l,...,n, (12) 
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in involution, {Fj, Fk} = 0. So, the system under consideration is completely 
integrable. 

The generating function of integrals of motion has the form 



IV IV £ IV 



2 n /2 

jk 



and we have 

{G x ,G^} = 0. 

Let us rewrite these equations in the Hamiltonian form. For this, we use 
the Poisson structure related to the Lie algebra e{n) of motion of Euclidean 
space M™, 

{J'iji l>km\ $ik Ijm $irn Ijk ^jk^irn ^jrn^iki 

(13) 

{hj, Vk} = $jk Vi ~ $ik Vj, {Uj, Uk} = 0. 

So, the quantities yj, lk m generate the Lie algebra e(n) of motion in n- 
dimensional Euclidean space. 

The Hamiltonian H should be the function of Fj and in fact, due to the 
quadratic form of equations of motion, it should be a linear combination of 

H = \b 3 F j . (14) 
The simple calculations show that bj = aJ 1 . Hence 
Vj = {H,Vj}, i jk = {H,l jk }, 



(15) 



Eqs. (15) are the Hamiltonian equations related to the Poisson structure for 
the Lie algebra e(n). 

For two-dimensional case these equations were integrated explicitly by 
Weber [5] in terms of theta functions of two variables. For multi-dimensional 
case, see [9] and [12]. 
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Note that as it was shown by Joachimsthal [13], the quantity 

I(x,y)=A(x)B(y) 

is an integral of motion, I(x, y) = c 2 . 

Now we may obtain the expressions for the quantities Xj and A(x): 



r) = - B-\y) a 3 A Vj , A(x) = ^ | = c 2 B~\y). (16) 



Finally, the expression for t as the function of r is given by a quadrature 

t( T ) = [ A{t) dr. (17) 
Jo 

So, we get the formulae for all dynamical variables. 
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Appendix 

Lettre de M. Jacobi a M. Arago, concernant les lignes geodesiques tracees 
sur un ellipsoide a trois axes [1]. 

Konigsberg, le 28 decembre 1838 
Monsieur, 

Je suis parvenu a ramener aux quadratures la ligne geodesique sur un 
ellipsoide a trois axes inegaux. Ce sont des formules extremement simples, des 
integrates abeliennes qui se changent dans les integrates elliptiques connues, 
en egalant entre eux deux de ces trois axes. Ce probleme m'ayant paru 
long-temps tres difficile, je crois que sa solution pourra interesser peut-etre 
quelques-uns des illustres membres de lAcademie des Sciences. 

Acknowledgements I am grateful to Max-Planck-Institut fur Mathematik 
in Bonn, where this paper was prepared, for the hospitality. 
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